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Abstract 



The dynamics of stability and collapse of a trapped atomic Bose-Einstein 
condensate (BEC) coupled to a molecular one is studied using the time- 
dependent Gross-Pitaevskii (GP) equation including a nonlinear interaction 
term which can transform two atoms into a molecule and vice versa. We find 
interesting oscillation of the number of atoms and molecules for a BEC of 
fixed mass. This oscillation is a consequence of continuous transformation in 
the condensate of two atoms into a molecule and vice versa. For the study 
of collapse an absorptive contact interaction and an imaginary quartic three- 
body recombination term are included in the GP equation. It is possible to 
have a collapse of one or both components when one or more of the nonlinear 
terms in the GP equation are attractive in nature, respectively. 

PACS Number(s): 03.75.Fi, 05.30.Jp 
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I. INTRODUCTION 



The experimental detection [0-0] of Bose-Einstein condensation (BEC) at ultralow tem- 
perature in dilute trapped alkali metal and hydrogen atoms has initiated intense theoretical 
activities in this subject 0-0. Many properties of the condensate are usually described 
by the nonlinear meanfield Gross-Pitaevskii (GP) equation [|T3| which can be considered to 
be a generalization of the Schrodinger equation in the presence of interatomic interaction 
leading to a nonlinear term. 

Two interesting features of BEC are (a) the collapse in the case of an attractive atomic 
interaction and (b) the possibility of the formation of a hybrid atomic and molecular 
condensate |6Hl0i coupled by coherent Raman transition |jl6| or Feshbach resonance [P,|l7 . 



Although there has been experimental observation of a coupled BEC involving two atomic 



states |]T8|, so far there has been no experimental realization of a molecular BEC. Because 
of the intrinsic interest in a hybrid atomic and molecular BEC, in this work the dynamics 
of such a BEC is studied using the time-dependent GP equation with special attention 
to its stability and collapse. For the case of a trapped two-species BEC a large repulsive 



interaction between the two species can induce instability in the system [19 



For attractive interatomic interaction [|^,^, the condensate is stable for a maximum 
critical number of atoms. When the number of atoms increases beyond the critical number, 
due to interatomic attraction the radius of BEC tends to zero and the central density tends 
to infinity. Consequently, the condensate collapses emitting atoms until the number of 
atoms is reduced below the critical number and a stable configuration is reached. With a 
supply of atoms from an external source the condensate can grow again and thus a series 
of collapses can take place, which was observed experimentally in the BEC of ^Li with 
attractive interaction [0,0. Theoretical analyses based on the GP equation also confirm the 
collapse 

The possibility of the formation of a hybrid trapped atomic BEC coupled to a BEC of 
diatomic molecules via either Raman photoassociation |T^ or a Feshbach resonance [^JT] 
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has been suggested and studied recently 0,01 • These authors used a nonhnear interaction 
term P,|20| in the Hamiltonian that allows for a coherent number-conserving nondissipative 



transition of two atoms to a molecule and vice versa. Possible candidates for molecular 
BEC are "^^Rbz, ^^Rba, ^^Na2 etc 0,|2T[. The coherently coupled BEC could allow localized 



solitons even in the absence of a trap [0. However, in this work we shall only consider the 
interesting dynamics of this hybrid BEC in the presence of a trap. Because of the general 
interest there have been several theoretical studies involving two coupled atomic condensates 
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Here we critically study the atomic plus molecular hybrid BEC in the presence of attrac- 
tive bosonic interaction(s). The usual CP equation allowing for transition of two atoms to a 
molecule and vice versa conserves the overall number of atoms. In the study of the dynamics 
of this BEC, oscillation of the individual numbers of atoms and molecules are found with the 
conservation of overall total number (number of atoms plus twice the number of diatomic 
molecules). The dynamics of collapse for attractive interaction is studied by introducing 
absorptive contact interactions which allow for a growth in the particle numbers from an 
external source in addition to imaginary quartic three-body interactions responsible for re- 
combination loss [Q. In the presence of the imaginary contact and three-body interactions 
the GP equation does not conserve the overall number of atoms. If the strengths (and signs) 
of these two imaginary interactions are properly chosen, the solution of the GP equation 
could produce a growth of the condensate with time for numbers of atoms and molecules 
less than the respective critical numbers. Once they increase past these critical numbers 
the three-body interaction takes control and these numbers suddenly drop below the critical 
level by recombination loss signaling a collapse 0]. Then the absorptive term takes over 
and the individual numbers start to increase. This continues indefinitely showing an infinite 
sequence of collapses. 

There are three types of nonlinear interactions in the GP equation for the hybrid BEC 
of atoms and molecules: e.g., between two atoms, between two molecules, and between 
an atom and a molecule. If one of these is attractive, one component of the condensate 



may experience collapse. If two are attractive one can have collapse in both components. 
Specifically, one can also have collapse of both components if only the interaction between 
an atom and a molecule is sufficiently attractive. 

In Sec. II we present a set of coupled GP equations which is appropriate for the study of 
the hybrid atomic-molecular BEG including a new nonlinear interaction term which allows 
for the possibility of transformation of two atoms in the condensate to a molecule and vice 
versa. In addition we include an absorptive contact term and a three-body recombination 
term appropriate for studying collapse of the BEG. In Sec. Ill we solve the coupled GP 
equation numerically and study the hybrid atomic-molecular BEG under different possibil- 
ities of atomic and molecular interactions. Finally, in Sec. IV we give some concluding 
remarks. 



II. COUPLED GROSS-PITAEVSKII EQUATION FOR ATOMIC-MOLECULAR 

CONDENSATE 

Here, in addition to the kinetic energies, trap potentials, and usual nonlinear interaction 
terms the following interaction responsible for a nondissipative transfer of two atoms into a 
molecule and vice versa is also considered in the hamiltonian |]6|,|20[ 



^int = f ' 



^i^i + ^r^2i , (1) 



where ipi {il)2) is the complex atomic (diatomic molecular) field and x the coupling. For a 
nonzero x neither the number of atoms nor that of the molecules is conserved. However the 
total mass of all the atoms and molecules is conserved. In a BEG the chemical conversion of 
two atoms to a diatomic molecule is dominated by coherent stimulated emission, in which 
transitions are enhanced by the number of molecules already occupying the ground state. 
The conversion of atoms to molecules could be resonant via a Feshbach resonance. Or, it 
could be effected via stimulated Raman coupling introduced by two laser fields. In this 
fashion it is possible to select a specific molecular state in the condensate. 
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In the presence of the extra interaction (|l]) the mean-field GP equation for wave-function 
components become 

2 



h"^ 1 8^ 1 

2 2 X — ^ I / / \i2 



d_ 



i^iir, r) + xi^2{r, T)^l{r, r) = 0, 



(2) 
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r + cmuj^r'^ + ^ g2j\ijj{r, r) 



V'2(r,r) + |^/'i(r,r) 



0, 



(3) 



where m is the atomic mass, 2m the diatomic molecular mass, Qij the nonlinear interaction, 
uj the frequency of the harmonic oscillator trap, and i = 1 (2) denotes an atom (molecule). 
The parameter c has been introduced in (|^) to modify the frequency of the trap for molecules. 
The diagonal interaction term gu is related to the scattering length of i-type bosons via 
gu = Anh^ai/mi, where m, is the corresponding mass. The wave-function components are 
normalized as 



d^r\Mr,r)\' = N, 
d'r\^2{r,T)\^ = N2 



(4) 
(5) 



where Nj is the number of j-type condensed boson. For x = the number of bosons of 
each type is conserved. For x 7^ 0, the number of bosons of each type is not separately 
conserved and there will be continuous transfer of one type of bosons to another and vice 
versa. However, the overall total number (A^^i + 2N2) is conserved. 

As in Refs. []13[ it is convenient to use dimensionless variables defined by x = y/^r/a^^^, 



t = Tuj, and (f)i{x,t) = xtlJi{r,T){^/27ra'^^y^^ where aj^Q = \J fi / [muj) . In terms of these, (0) 
and (13) become 



^x'J^ |0,(x,t)p \Mx,t)r 

d 



dx 



x^ 



x" 



+ 171-1 



dt 



01 (a;, t) + T] 



(t)2{x,t)<p\{x,t) 



0, 



X 



(6) 



02(a;,t) + - — — = 0, (7) 



where riij = gijm/ {y^Trfi^a-^Q) is negative (positive) when the corresponding interaction is 
attractive (repulsive) and the dimensionless couphng i] = x/[/iC(j(-\/2vra|^^)^/2], xhe diagonal 
elements of Uij are ria = 2v^(aj/aj^Q). In (|^) and (|^) a diagonal absorptive (strength 7^) and 
quartic three-body (strength ^j) interactions, appropriate to study collapse, are introduced 
A nonzero 7^ allows the possibility of the absorption of atoms and molecules from the 
surrounding into the condensate, whereas a nonzero allows for the possibility of ejection 
of particles from the condensate due to three-body recombination. Consequently, the total 
mass of the condensate is not conserved for nonzero 7^ and A nonzero 7^ and C,i is 
appropriate for the study of collapse of the system. 

The fluctuation of atomic and molecular numbers is best studied via the quantities 

1 

^at = -(5)/ \M^,t)\'dx, (8) 
1 

^mol = TToT / \Mx,t)\'^dx, (9) 



ATf) Jo 

where Ni^^ and A'"2°^ are the initial numbers of the two types. In the presence of a general 
absorptive interaction 77 7^ 0,7^ 7^ and C,i 7^ 0, the quantities N^^ and N^^^ carry the 
information about time-evolution of the number of atoms and molecules in the condensate. 
For ji = C,i = 1] = 0, the number of atoms and molecules are separately conserved and one 
has Ng^^ = A^jj^qJ = 1. For 7i = = and ?7 7^ 0, there is no separate number conservation 
for atoms and molecules: N^^^ 7^ 1, N^^^ 7^ 1; although the total mass of the condensate 
is conserved. For 7^ 7^ and/or C,i 7^ there is no conservation of the total mass of the 
condensate. 

The root-mean-square (rms) radius Xrms(^) of component i is given by 



3;rms(^) 



j,^ x'Mx,t)\V^' 



J,^Mx,t)\^dx 

The oscillation of this rms radius denote radial oscillation of the condensate. 
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(10) 



III. NUMERICAL METHOD 



To solve and (|^) numerically one needs the proper boundary conditions as x ^ and 
oo. For a confined condensate, for a sufficiently large x, (f)i{x,t) must vanish asymptotically. 
Hence all the nonlinear terms proportional to the product of wave-function components can 
eventually be neglected in the GP equation for large x. Consequently the the asymptotic 
form of the physically acceptable solution is given by lima;^oo \4'i{x,t) \ ~ exp(— Next 
we consider and (0) as a; 0. The nonlinear terms approach a constant in this limit 
because of the regularity of the wave function at x = and one has the condition |0i(O, t) \ = 
0. 

Next we discretize (||) and (0) in both space and time by using a space step h and time 



step A with a finite difference scheme using the unknown 



which will be approximation 



of the exact solution 0j(a;p,tfc) where Xp = ph and = kA. After discretization we obtain 
a set of algebraic equations which could then be solved by using the known asymptotic 
boundary conditions. The numerical procedure is similar to that in refs. The time 

derivative in these equations involves the wave functions at times tk and t^+i- As in the 
uncoupled case we express the wave functions and their derivatives in (j^) and (|^) in terms 



of the average over times and t^+i ||2^ and the resultant scheme leads to accurate results 
and good convergence. In practice we use the following Crank-Nicholson-type scheme |23 
to discretize the partial differential equations (||) and 
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Considering that the wave function components (pi are known at time t^, (|TTD and (|T2D 
are equations in the unknowns — (0i)p+i, i(f)i)p'^^ and (0i)pli. In a lattice of points (|ll]) 
and ( p!2|) represent a tridiagonal set for p = 2, 3, (A^— 1) for a specific component 0j. This 



set has a unique solution if the wave functions at the two end points and (0i)^^ 



are known |£3|. In the present problem these values at the end points are provided by 
the known asymptotic boundary conditions. This tridiagonal set of equations is solved by 



Gaussian elimination and back-substitution 23 



To solve (pJ]) and (0) we employ space step h = 0.0001 with xmax < 15 and time step 
A =0.05 with tmax = 100. After some experimentation we find that these values of the 
steps give good convergence. The iteration is started with the known normalized (harmonic 
oscillator) solution of @ and (0) obtained with riij = 'Ji = = rj = with normalization 

r \Mx,t)\' = Ni'\ z = l,2, (13) 



where A'^''^^ is the initial number of atoms. The nonlinear parameters riij and rj are increased 
by equal amounts over 500 to 1000 time iterations starting from zero until the desired final 
values are reached. During this process we keep = 7i = 0. The resulting solution is the 
ground state of the condensate corresponding to the specific nonlinearity. 

The time- dependent approach is the most suitable for solving time-evolution problems. 
In the present study we only consider evolution problems starting from a stable condensate 
at t = 0. In these cases the stationary problem is solved first and the wave function so 
obtained serves as the starting wave function for the time-evolution problem. 

To study the effect of an external source of bosons on preformed condensate (s) we need 
to set 7j 7^ and ^ 7^ and perform the time iteration of the GP equations (|^) and (|^ 
starting at t = from the solution obtained with 7^ = and C,i = 0. A nonzero 7^ will 
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make the condensate grow in size absorbing bosons from outside. In the case of attractive 
interaction(s), as the condensate grows past the critical size a nonzero ^ will allow for three- 
body recombination so that the condensate can shrink and experience the so called collapse 
by emitting particles. Hence a nonzero 7^ and could simulate the dynamical growth of 
the condensate in general as well as its collapse for attractive interaction. 

IV. NUMERICAL RESULT 

Instead of concentrating on a specific atomic-molecular BEC, here we study some uni- 
versal features of such a BEC related to stability and collapse. For a specific application we 
require the knowledge of the different scattering lengths: those for two atoms, two molecules 
and between an atom and a molecule. At present these scattering lengths are not all known 
for a specific system. 

A. Repulsive Interaction: Stable Condensate 

First we consider the nondissipative solutions of (^) and obtained with r/ 7^ and 
7j = = for all repulsive interactions and demonstrate the oscillation of numbers of 
atoms and molecules. In this case the condensate increases in size indefinitely with time 
if a nonzero 7^ and are used. This is why it can be termed stable. However, there is 
spontaneous transformation of atoms to molecules and vice versa even for 7i = = and 
the number of atoms and molecules oscillate with time. The numerical results are convergent 
for small values of nonlinear parameters riij and 77. Here to demonstrate the viability of the 
present numerical scheme we present results for reasonably large values of n^. Specifically, we 
consider c = 1/16, ivf^ = A^f^ = 10000, riu = n<n = 0.001, = n^x = 0.0005, Xmax = 15 
and 1] = 0.001, 0.005, and 0.01. For a typical experimental situation ||11| this corresponds 
to 2^/2\ai\ / a-^Q = 10'^. The solution of the complete coupled GP equations (||) and (|^ 
is obtained by iteration starting with the solution of the uncoupled Schrodinger equation 
for the harmonic oscillator problem obtained by setting Uij = 77 = in these equations. 
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The nonlinear parameters riij are increased in small steps during iteration until their final 
values are reached after 1000 iterations. We then obtain the desired solution of the coupled 
GP equation. In figure 1 (a) we plot the atomic and molecular wave functions 4>i and 4>2 
obtained from this procedure. 

We continue the time iteration after the solutions of the coupled GP equation shown in 
figure 1 (a) are obtained at a specific time, which we now call t = 0, without changing the 
nonlinear parameters anymore. In figure 1 (b) we plot the quantities N^^ and -?Vjjiol 
time t ior t > 0. For 77 7^ 0, as shown in figure 1 (b), the numbers Nq^^ and N^^^ undergo 
rapid oscillations. In figure 1 (c) we plot the corresponding rms radii of the atomic and 
molecular condensates for f > 0. We find that the rms radii also execute oscillation along 
with A'^g^^ and N^^^ due to continued dynamical conversion of atoms to molecules and vice 
versa. The oscillation of the rms radii is related to the radial oscillation of the condensate. 

B. Attractive Interaction: Route to Collapse 

Now we consider the case of attractive atomic and molecular interaction (s) which may 
lead to collapse. We note that is a nonlinear parameter in the coupled GP equation and a 
large value of the nonlinear parameter 77 will imply a large numerical error in the solution. A 
very small rj on the other hand will not give adequate feedback between atoms and molecules. 
We study collapse in the hybrid BEG under several conditions with r] = 0.00625, as this 
value of T] gives adequate feedback between atoms and molecules and accurate numerical 
result. 

First we study the simplest case of collapse by taking all interactions to be attractive 
corresponding to a negative riij for all i and j. Now both components of BEG could experi- 
ence collapse. In this case the wave function components calculated numerically as in figure 
1 (a) at t = are plotted in figure 2 (a) for A^f) = = 1600, nn = -0.000975, 7122 = 
-0.000475,7112 = 7121 = -0.000125, c = 1/16, 77 = 0.00625, xmax = 15, with = 7i = 0. 
In the case of attractive interaction, for BEG with ^Li the typical number of atoms in the 
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condensate is 10^ and |ai|/%Q ~ 10~^ to 10~^ In this example, |ai|/a]^Q ~ 3.5 x 10~^, 
and |a2|/ctiio ~ ^'"^ ^ 10^^. Hence the present numbers are well within the experimental 
scenario of BEC for atoms with attractive interaction. From figure 2 (a) we find that the 
wave-function components are more centrally peaked compared to those in figure 1 (a). 
This corresponds to small rms radii and large central density denoting an approximation to 
collapse. 

Although the collapse of the coupled condensates could be inferred from the wave func- 
tions of figure 2 (a), we also study the dynamics of the collapse from a time evolution of 
GP equations @ and (0) for t > in the presence of an absorption and three-body recom- 
bination, e.g., for 7i 7^ and 7^ @]. We consider the solution of (D and as shown 
in figure 2 (a) at t = and allow this solution to evolve in time with 7i 7^ and 7^ 
by iteration. The general nature of time evolution is independent of the actual values of 
7i and C,i employed, provided that a very small value for 10~^) and a relatively larger 
one for 7j(~ 0.1) are chosen [^]. In actual calculation we took 71 = 0.15,72 = 0.08 and 
^1 = 0.00000031,^2 = 0.00000094. The quantities A^^t and N^^^ show the succession of 
collapses in a transparent fashion and we plot them in figure 2 (b) vs. time t. These quan- 
tities increase with time due to the absorptive interaction (74). Once they are larger than 
the respective critical numbers for collapse, the three-body recombination takes over and 
these quantities fall below the critical level and the system collapses. Then the absorptive 
term is operative and these quantities increase again resulting in a sequence of collapses. 
In figure 2 (c) we plot the rms radii Xrms of the atomic and molecular condensates of the 
model of figure 2 (b) undergoing collapse for t > 0. The rms radii of the components execute 
oscillation similar to those in the uncoupled case 

The collapse was studied under various situations. After numerical investigation we find 
if only the diagonal parameter rin or 77-22 is attractive and all other repulsive then it 

is possible to have a collapse of component 1 or 2. If in addition to one diagonal parameter 
riii, the nondiagonal parameters ni2 and n2i are attractive, it is possible to have a collapse 
in both components. If both diagonal parameters nu and n22 are attractive, one can have 
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collapse of both the components. 

Instead of reporting all our studies, in the following we consider the special case of collapse 
in the case where the diagonal elements rin and n22 are positive and the nondiagonal terms 
ni2 and ^21 are negative. This corresponds to repulsion between two atoms and between two 
molecules and an attraction between a molecule and an atom. This case is of relevance, as 
experimentally, the rubidium system has been measured to have strongly attractive atom- 



molecule interaction and repulsive atom-atom and molecule-molecule interactions P,|3,|T0 
Also, it is now possible to vary continuously the nature (attractive or repulsive) and strength 
of interaction between bosons p4| which can create a specific situation of atom- molecule, 
atom-atom and molecule-molecule interactions. We find that when ni2 and 7221 are suffi- 
ciently attractive, they can surplus the repulsive terms nn and n22 in and (|^ and the 
overall effective contribution of the quadratic nonlinear terms in these equations becomes 
attractive and one can have a collapse of both the components. To illustrate this interesting 
effect we study numerically the collapse for a negative ni2 and n2i and positive rin and 
^22. In figure 3 (a) we plot the wave-function components obtained by solving and (|^ 
with n[°'^ = = 1600, nu = 0.0005,^22 = 0.000625, = ^21 = -0.00144, c = 1/16, 
T] = 0.00625, Xmax = 10 and 7i = 6 = 0. Both wave-function components are peaked near 
X = and have small rms radii. The system collapses with a small increase in |ni2| and/or 
|?T,2i|. The time evolution of Ng^^ and N^^^^^ and the corresponding rms radii are found by 
solving d^) and iteratively with nonzero C,i and 7, starting with the wave function pre- 
sented in figure 3 (a) at t = 0. For t > we employ ^1 = 0.00000031,^2 = 0.00000075 and 
7i = 0.12,72 = 0.05. The quantities N^^^ and N^^^ show the dynamics of collapse directly 
and they are plotted in figure 3 (b). The successive growth and decay of these numbers with 
time exhibit the sequence of collapses. Finally, in figure 3 (c) we plot the evolution of the 
corresponding rms radii. The results for < t < 100 in figures 2 and 3 are calculated with 
2000 iterations of the GP equation (|^) and (|^ using a time step 0.05. 



12 



V. CONCLUSION 



To conclude, we studied the stability and collapse of a trapped hybrid BEC consisting 
of atoms and molecules |]6|,[7|,PJT0[1 coupled by coherent Raman transition [16| or Feshbach 
resonance HJT7[] using the time-dependent coupled GP equation. An additional nonlinear 
interaction between the atoms and molecules given by (|^) is introduced in the formulation 
which allows spontaneous transformation of two atoms into a molecule and vice versa. The 
effect of an external source of boson on the condensate is simulated by the absorptive term 
7j and the effect of three-body recombination necessary for the study of collapse is included 
via the quartic nonlinear term C,i in the GP equation. 

We find interesting oscillation of the number of atoms and molecules due to the trans- 
formation of atoms to molecules and vice versa even when all interactions are repulsive. 
The component i of the condensate could experience collapse when the interaction among 
bosons of type i is attractive. Both components could experience collapse when at least the 
interaction between a molecule and an atom is attractive. They could also experience col- 
lapse when two of the bosonic interactions are attractive. The time evolution of collapse is 
studied via the time- dependent GP equation with absorption and three-body recombination. 
The number of particles of the component (s) of BEC experiencing collapse alternately grows 
and decays with time. The rms radii of the collapsing condensate also exhibit interesting 
oscillation as in the uncoupled case 0]. The oscillation of the rms radii denotes a radial 
oscillation of the condensates which can be observed experimentally. Recently, we studied 
the possibility of collapse in a coupled atomic BEC without any interaction allowing for 



transformation of one type of atom to another ||25|. With the possibility of observation of 



coupled BEC, the findings of this study and of ref. [|25[ could be verified experimentally in 
the future. 
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Figure Caption: 

1. (a) Atomic and molecular wave function components (l>i{x) and 02 (a^) vs. x at t — 0, 
(b) the corresponding numbers N^^ and -/Vj^q^ vs. time t {> 0), and (c) the rms radii vs. 
time t (> 0) obtained for A^f) = N^^^ = 10000, rin = 7122 = 0.001, ni2 = 7121 = 0.0005, 
c— 1/16, for 7j 0. The curves are labelled by the corresponding values of 77. 

2. The same as in figure 1 with A^f) = N^^'^ = 1600, nn = -0.000975, n22 = -0.000475, 
= 7121 = -0.000125, c = 1/16, and 77 = 0.00625. In (a) = 7i = 0, and in (b) and (c) 

^1 = 0.00000031, 71 = 0.15, 6 = 0.00000094, 72 = 0.08. 

3. The same as in figure 1 with A^f^ = A^f^ = 1600, 7iii = 0.0005, 7122 = 0.000625, 
= 7i2i = -0.00144, c = 1/16, and 77 = 0.00625. In (a) = 7^ = 0, and in (b) and (c) 

^1 = 0.00000031, 7i = 0.12, ^2 = 0.00000075, 72 = 0.05. 
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